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We derive an effective Hamiltonian for the ionic Hubbard model at half filling, extended to include nearest- 
neighbor repulsion. Using a spin-particle transformation, the effective model is mapped onto simple spin-1 
models in two particular cases. Using another spin-particle transformation, a slightly modified model is mapped 
into an SU(3) antiferromagnetic Heisenberg model whose exact ground state is known to be spontaneously 
dimerized. From the effective models several properties of the dimerized phase are discussed, like ferroelectric- 
ity and fractional charge excitations. Using bosonization and recent developments in the theory of macroscopic 
polarization, we show that the polarization is proportional to the charge of the elementary excitations. 
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I. INTRODUCTION 

The ionic Hubbard model (IHM) has been proposed 01 0] 
for the description of the neutral-ionic transition in mixed- 
stack charge-transfer organic crystals tetrathiafulvalene-p- 
chloranil.|301 In the 90's the interest on the model increased 
due to its potential application to ferroelectric perovskites.0 
IS 13 S HM E3 El This model is a Hubbard Hamiltonian 
with nearest-neighbor hopping t and on-site repulsion U, on a 
bipartite lattice with alternating diagonal energy ±| A for the 
two sublattices. At half filling and in the atomic limit it — > 0), 
the ground state is a band insulator (BI) for U < A (the sites 
with diagonal energy — - A are doubly occupied), but is a Mott 
insulator (MI) for U > A (all sites singly occupied). The sit- 
uation is similar to the extended Hubbard model (EHM) with 
A = but nearest-neighbor repulsion V. fl2l fo fl4ll In the 
EHM for t = 0, as U increases, there is a transition from a 
charge density wave (CDW) insulator, with alternating dou- 
bly occupied and empty sites to a MI for U = 2V. However, 
while in this model, the result persists for finite small t as it 
has been shown in second order perturbation theory in t. f\3\ 
perturbation theory in the IHM becomes invalid at U = A and 
non-trivial charge fluctuations persist even in the strong cou- 
pling limit. Furthermore while previous results in the IHM 
were interpreted in terms of only one transition, Fabrizio et 
al. used arguments based on field theory, which assumes weak 
coupling (essentially A <C U <C t), to propose an alternative 
scenario that includes two transitions as U is increased: the 
first one is an Ising-like charge transition at U = U c from the 
BI to a bond-ordered insulator (BOI). When U is further in- 
creased, the spin gap vanishes at U s > U c due to a Kosterlitz- 
Thouless transition between the BOI and the MI llOUllll An 
intermediate BOI phase also takes place in the EHM for large 
enough t. til 111 

After the proposal of Fabrizio et al. , several studies of the 
ID IHM tried to identify the number and the nature of the 
different phases. Torio et al. determined the phase diagram 
using the method of crossing of energy levels, which for this 
model turns out to be equivalent to the method of jumps of 
Berry phases. 11211 The basic idea is that while in finite sys- 



tems conventional order parameters such as charge and spin 
structure factors vary continuously at a phase transition, in 
a system of length L one may define charge and spin topo- 
logical numbers which change discontinuously at U C (L) and 
U S (L). The charge (spin) Berry phase has a step at U C (L) 
(U S (L)), and extrapolating these numbers to the thermody- 
namic limits leads usually to accurate results for the transition 
lines. El El u c (l) might be detected in transport experi- 
ments through some molecules or nanodevices. Jl8ll Torio et 
al. found a phase diagram consistent with the proposal of 
Fabrizio et al. for weak coupling. They also found that for 
strong coupling the intermediate BOI phase persists with a 
width U s — U c = 0.6£. Recently, similar methods were used 
to check the characteristics of each quantum phase transition 
and confirm previous results in systems up to 18 sites. B19I1 
However, Monte Carlo calculations for the bond-order param- 
eter, polarization and localization length were interpreted in 
terms of only one transition at U c and two phases (BI and 
In addition, some density-matrix renormalization- 
group (DMRG) calculations have not detected the transition 
at Z7 . ll2lll Instead, other DMRG calculations using careful 
finite- size scaling are consistent with three phases and two 
transitions and with the phase diagram of Torio et al. , although 
they suggest a smaller width of the BOI phase for small t 
(U s - U c ^ O.4t).0E3 

The numerical difficulties for detecting the Kosterlitz- 
Thouless transition are originated in the exponential closing 
of the spin gap as U approaches to U s from below. Conse- 
quently, direct numerical calculations of the relevant corre- 
lation functions are unable to detect a sharp transition at U s 
unless a very careful finite-size scaling is made. The same dif- 
ficulty appears in the Hubbard chain with correlated hopping, 
where the existence and the position of the transition is con- 
firmed by field theory iflfS l24ll and exact ll25ll results. In the 
scenario with only one transition there is still a question about 
the nature of the large U phase (U > U c ). Wilkens and Mar- 
tin 112 (ll suggested that the MI phase does not exist for finite t, 
i.e., the spin gap remains open for any finite U. This is in con- 
tradiction with the strong coupling expansion for t <C U — A, 
which maps the IHM onto an effective spin Hamiltonian H 
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with closed spin gap. 

nuns 

The absence of a spin gap is 
also supported by the power-law decay of charge-charge cor- 
relation functions (in spite of the presence of a charge gap) 
that have been calculated numerically and analytically 
(using H). $2(h with very good agreement between both re- 
sults. This is a consequence of the strong charge- spin coupling 
in the low-energy effective theory. These results suggest that 
if there is only one transition, the BOI phase should disappear 
from the phase diagram. We mention here that in presence of 
electron-phonon interaction, the width of the BOI phase in- 
creases and, in the adiabatic approximation, it replaces the MI 
phase. GSeSi 

From the above discussion, it is clear that the existence 
and extension of the BOI phase in the electronic model de- 
serves further study, particularly in the strong coupling limit, 
t <C [/, A, since the arguments of Fabrizio et ah are not di- 
rectly applicable. In addition, little attention has been given 
so far to the electric polarization or ferroelectric properties 
of the BOI phase. |23 1231 The BOI phase is an electroni- 
cally induced Peierls instability that generates a ferroelectric 
state. Experimentally, a bond ordered ferroelectric state was 
observed in the pressure-temperature phase diagram of the 
prototype compound, tetrathiafulvalene-p-chloranilj3l 01 . In 
addition, as it was pointed out by Egami et al, |5|] the mi- 
croscopic origin of the displacive-type ferroelectric transition 
in covalent perovskite oxides like BaTi03 is still unclear. It 
has been known for a long time that a picture based on static 
Coulomb interactions and the simple shell model is inade- 
quate to describe some ferroelectric properties. O In a re- 
cent paper. ll27ll using a spin-particle transformation, 
we mapped a Hamiltonian which is close to the one obtained 
in the strong coupling limit of the IHM (in a sense that will 
become clear in Section IV) into an anisotropic SU(3) anti- 
ferromagnetic Heisenberg model that becomes isotropic and 
exactly solvable in ID for certain combinations of the pa- 
rameters. The exact solution is a spontaneously dimerized 
BOI that helps to understand the physics of the IHM in the 
strong coupling limit. The large value of the correlation length 
£ = 21.0728505... allows to understand the numerical diffi- 
culties for detecting this phase in finite size systems. We also 
used a second spin-particle transformation B31I1 that maps the 
constrained fermions into 5 = 1 SU(2) spins and brings the 
effective model into a biquadratic Heisenberg model with a 
spatial anisotropy along the z-axis which is proportional to 
A — U. This allows to interpret both transitions and the frac- 
tional charge excitations of the BOI phase in the simple lan- 
guage of S = 1 spins. 

In this work, we derive an effective Hamiltonian H e f f for 
the strong coupling limit of the extended IHM (including A 
and V). As in the t — J model, the single-site Hilbert space of 
H e ff has three states. In two particular relevant cases H e ff 
is mapped into simple models for S = 1 SU(2) spins. This 
opens the possibility of studying both quantum phase tran- 
sitions in ID with the quantum loop Monte Carlo technique 
which is specially designed to deal with criticality due to the 
non-local nature of its dynamics. 1 32] Harada and Kawashima 
ll33ll showed that this technique provides an accurate determi- 
nation of the critical parameters of a Kosterlitz-Thouless tran- 
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sition. The mapping to the exact solution l27ll is briefly re- 
viewed. We discuss the polarization of the BOI phase. In par- 
ticular, we construct the bosonized expressions for the "twist" 
operators z c L and z s L 10, |3 |35|l which gives information on 
the conducting and polarization properties of the system, and 
are also related with the Berry phases. This allows us to re- 
late the polarization in the BOI phase with the charge of its 
elementary excitations. 

II. THE EFFECTIVE HAMILTONIAN 

In this Section, we derive an effective Hamiltonian for the 
extended IHM (EIHM) valid for t <C A, but any value of U 
if V is small. 13611 For simplicity we restrict our study to ID. 
The extension of this derivation to dimension higher than one 
is straightforward but somewhat cumbersome for The 
Hamiltonian of the EIHM can be written in the form 

H = -iB/U^ + H-cJ + ^Q-l)^ 

i,<j i 

+ U ^(n iT - i)(n a - 

i 

+V^(n;-l)(n; +1 -l), (1) 

i 

where n iG = f- a f ia , ra* = J2 F n ia> and A,£ > 0. This 
Hamiltonian has four states per site: empty, doubly occupied, 
or singly occupied with spin up or down. At half filling and 
for t = 0, all odd sites (those with energy — A/2) have at 
least one particle, and no even site is doubly occupied. Thus, 
the relevant low-energy Hilbert subspace Ho has three states 
per site. Our aim is to eliminate linear terms in t which mix 
states of Ho with the rest of the Hilbert space. This can be 
done with a canonical transformation that leads to an effective 
Hamiltonian H e f f which acts on Ho and contains terms up to 
second order in t. The procedure is completely standard (see 
for instance Ref. To simplify the construction of H e f / , 

it is convenient to perform an electron-hole transformation for 
the odd sites 

4 T = ~fil> 4[ =fiV foroddz 

4a = fli for even z. (2) 

The Hamiltonian now becomes invariant under translation of 
one site (i — > i ± 1) but it does not conserve the number of 
particles due to the particle-hole transformation. The U(l) 
symmetry associated with the conservation of the total num- 
ber of electrons is now generated by a the staggered charge 
operator Q = J2i a (~ ^T c la c ia- The transformed Hamilto- 
nian becomes 

i,a i 

-V^-ljf^-lJ + f/^n^. (3) 

i i 
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Now, for t <C A but arbitrary A — U. l^k we can eliminate 
the states with double occupancy on any site. This constraint 
can be incorporated to the fermionic algebra by defining the 
constrained fermion operators: 



4a 



rt_ c - 



), 



■ r )- (4) 



Performing a canonical transformation we eliminate the terms 
proportional to t. These terms mix the low energy subspace 
Ho with the orthogonal high-energy subspace Hi, which con- 
sists of all the states having at least one doubly occupied site 
in this representation. Up to second order in t, we obtain the 
effective Hamiltonian: 



T C i+l| 



'I+1T 



H.c.) 



V ^(1 - m)(l - rn+i) + E ^(n, - 1) 



SJ = V2(K] 
S] = n jT - n ji 



3 C 3\ 



+ 4 



'3 1 



where K, is the kink operator l3lll : 



exp[z7r 



k<j 



(7) 



(8) 



One way to visualize the effects of this transformation is to 
note the corresponding mapping of the fermion states at any 
site to the spin states \SS Z ): \j0) -> |jl0), c^Kj\0) -> 

|jll), and Cj^Kj\0) — > \jl — 1). Up to an irrelevant con- 
stant, the resulting i^ e j / to first order in t is 



+ ^ J^Sj • Sj+i - -njUj+i) 

+ E ^l+^-<^(y + 2s,-s,_,), (5) 

where E = (A — U)/2, and is the spin vector at site z, the 
components of which are sf = 1/2 ^ r r , cJ T cr" r ,c ir / with 
a = {x,y,z} {o v are the Pauli matrices). The exchange in- 
teraction Ji and correlated hopping if 1 coefficients are oper- 
ators if V ^ 0. If V = 0, Ji = 2t 2 /(U + A) and tf = t 2 / A 
are constants. Ji comes from a second order process in which 
in the original Hamiltonian, the virtual intermediate state has 
a doubly occupied site at energy -f A/2 and an empty site at 
energy — A/2. In the original language, if 1 represent a cor- 
related hopping between two second nearest-neighbor sites at 
energy + A/2(— A/2) involving an intermediate state with an 
empty (doubly occupied site) at energy — A/2(+A/2). 

For general V, the expressions of these operators are the 
following: 



Ji 2t 



ni-i + n i+2 ~ 2ni_in i+2 



U + A-V 
(1 - nj_i)(l - n i+2 ) 
U+A+V 1 



t 2 l-rni_ 2 ni+2 



U + A 



(6) 



V 



(1 - ^_ 2 )(1 - n i+2 ) 



+ + A ^ y )(ni-2 + n i+2 ~ 2n i _ 2 n i+2 )/2]. 

While the general expression of H e f f is rather complicated, 
it takes simpler forms for some particular cases. For instance, 
if we only keep the linear terms in t, what amounts to neglect 
Ji and if 1 , the resulting already contains the non-trivial 
charge fluctuations of the EIHM near both transitions and is 
the minimal model to describe the system in the strong cou- 
pling limit. We can rewrite H£f f as a spin 1 model by using 
the spin-particle transformation 13111 



eff 



E 



[t(sfsf +1 + sf sy +1 )Pi^i 

V{Sff{Sl +1 f + (E + 2V)(S?)% (9) 



where Pij projects over the states with Sf + *S| = (one may 
use = 1 - (Sf + S*) 2 in Eq. Q ). 

When the second order terms are included for V = (as in 
the original IHM), the effective spin Hamiltonian becomes: 



H eff = T,[t{SfSf +1 + SfSf +1 )P iti+1 + E{Sff 
{{(SfS- +l ) 2 + {Sf +1 S-) 2 + 2Sf S? +1 - 2(Sf ) 2 (5f +1 ) 2 } 

^ch \ r ( qz \2 fox qx | qV qU \{qz \2 



S=±l 



x{(StS-_ 5 ) 2 + (S+_ 5 S-) 2 

+2SfSf_ 5 + 2(Sf) 2 (Sf_ 5 ) 2 }], 



(10) 



where J = 2t 2 /(^" + A) and t ch = t 2 /A. 

i^gjj is the minimal model that describes the physics of 
the EIHM in the strong coupling limit. It can be studied with 
quantum Monte Carlo or DMRG and its solution can bring 
further insight into the nature of the BOI phase, the precise 
position of the phase transitions, and the effect of V on them. 
Note that these spin 1 Hamiltonians as well as Eq. (l27l be- 
low, conserve not only total S z , but also the original charge 
Q = £) i (-l) < (S , ?) 2 |, and the ground state is in the Q = 
subspace. This seems to be an important reduction of the rel- 
evant Hilbert space for numerical calculations. 



III. LIMIT OF LARGE V 

For V ^> t, the ground state of the EIHM can be ob- 
tained by second-order perturbation theory in t, extending pre- 
vious results for the EHM. tT^I One can start either from H 
or H e ff. In the latter case, one can easily see that if 1 be- 
comes irrelevant and that there is an additional contribution 
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J a d = 2t 2 /(U — A — V) to J when empty sites are elimi- 
nated by another canonical transformation. 

Starting from a CDW with maximum order parameter, the 
energy per site up to second order in t becomes 



-E-V 



t 2 



(11) 



where the subscript II means "ionic insulator" to distinguish 
it from the BI because correlation effects due to V are also 
present. In the MI phase, one has an effective Heisenberg 
model with exchange interaction J tot given by Eq.@ for all 
ni = 1 plus J a d- From the exact solution, the ground state 
energy of this model is — Jtot In 2, and therefore 



where c (a) denotes charge (spin), H v (y = c, a) is a sine- 
Gordon Hamiltonian 



H " = J dx{^[7rK u Ul(x) 



ttK v 



+ cos(\/8</>„)}, 



and the charge-spin interaction is 



(14) 



/2A 
dx — cosfv^c) cos(v // 20 cr ). (15) 
na 

In terms of the density operators v kr for particles with wave 
vector k moving in the r (+ or -) direction, the phase fields are 




E 



MI 



-Jtot In 2; J+. 



4At 2 



[(U - V) 2 - A 2 ] 



(12) 



For t = 0, the transition between both phases is at U = 
A + 2 V. At this point, when V ^> t, the perturbative correc- 
tions in t 2 converge and the charge fluctuations of both phases 
have an energy cost V ^> t. Consequently, the essential in- 
gredient which can eventually lead to the formation of the BOI 
phase is absent. Therefore, there are only two phases and one 
transition in this limit. 



IV. BOSONIZATION 

The previous sections and most of the paper are concerned 
with the strong-coupling limit, in which t is small com- 
pared either to A or the interactions U, V. When the in- 
teractions are small, usually the continuum limit field the- 
ory plus renormalization group allows a precise description 
of the system. ll38H39ll In the EIHM, the usual procedures fail 
because A is a strongly relevant operator. Nevertheless the 
bosonized Hamiltonian was used to infer the existence of the 
intermediate BOI phase in the IHM and the character of its 
low-energy excitations. Here we describe the theory 

for the EIHM plus exchange interaction to draw conclusions 
to be used later, and construct the bosonized version of the op- 



erators U F S = exp[i^ J2j J a ( n j] =1= n jl)\> which enter the 
quantities z c £ s = (g\U^ s \g) where \g) is the ground state, L 
the length of the system and a the lattice parameter (and short 
distance cutoff below). 



A. Hamiltonian and main results 

The bosonized form of H + J J2i s i ' s i+i is 



H = H c + H a + H ccn 



(13) 



(j) v {x) 
where 



N,7TX 



-r l \k\ X -ik X[vk+ + Vk _ ) 



J v 



(16) 



[(7V +T + 7V_ T ) db (7V +T + M T )]/>/2, 



J Cj a = [(JV+t - A_ T ) db (7V +T - A- T )]/V2, (17) 

and N ra is the operator of the total number of particles with 
spin a. 

The parameters entering Eqs. (fTH) are determined by the 
following equations in terms of the usual g-ology coupling 
constants: 



04|| + #4± 



2ta- 



9c 



9l\\ ~ 2 #2|| - 2^ 2 ±, 9a 



2ta- 
= 9i\ 



9M\ 



9i± 
2n 2 '' 

94± 



292 



\2nv v F ■ 



9v 



2ttv f 



9v 



2# 2 ±, 
(18) 



The masses entering Eqs. (IT4t are proportional to the Umk- 
lapp coupling g 3 ± = a(U — 2V + 3J/2) and the backward- 
scattering one gi± = g a = a(U — 2V — J/2), while g c < 
for positive U, V, and J. 11611 For A = 0, the charge and spin 
sectors are decoupled and the usual renormalization group 
procedure for the sine Gordon model can be applied. By in- 
creasing U, there is a transition from the II (CDW) to the BOI 
at U c = 2V - 3 J/2 and another one from the BOI to the MI 
at U s = 2V + J/2.\ 16] Including second-order corrections to 
the coupling constants, it has been shown that for J = the 
BOI phase still exists for weak U ~ 2^.(3 The sign of g%± 
(#i_l) determines the charge (spin) sector. The width of the 
BOI phase goes to zero with J. 
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For A ^ 0, V = J = 0, Fabrizio et al | 10] replaced H c by 
an effective free massive theory for the MI phase. Then, treat- 
ing H ca perturbatively, they obtained an effective H a with a 
renormalized m a , showing that lowering U, a spin gap opens 
before the charge transition and therefore a BOI phase exists 
between the two transitions. They also discussed the exci- 
tations of the BOI phase taking the interaction part of the 
Hamiltonian as a phenomenological Ginzburg-Landau energy 
functional F with effective interactions. For its importance in 
what follows we briefly review this analysis. Ignoring some 
constants and factors and dropping the subscript _L, we can 
write 

a v = cos(v / 2^). (19) 

In Eq. JT9l and in the rest of this subsection, the gi and A 
are effective parameters renormalized by the interactions and 
their numerical values are different from those given previ- 
ously. Since we are in the spin sector of the BOI phase, this 
implies that the spin gap is open and therefore the effective 
gi < 0. Then, \a s \ = 1 minimizes the energy. If also 
gs < (corresponding to the BI phase) the energy is mini- 
mized by either a c = 1, a s = — 1 or a c = — 1, a s = 1. If 
one has a soliton (topological excitation) in the system be- 
tween these two vacua, clearly the jump in both phases is 
A(j) v = ±7r/\/2. Since the smooth part of the charge p(x) 
(spin S z (x)) operator is — V2d x (j) c /7r (d x <fi s /(V27r)) this ex- 
citation has total charge 1 and spin 1/2, as expected. If 
gs > 2A, the minima occur for a c = A/(2gs), a s = —1 
or a c = —A/(2gs), a s = 1. For an excitation between vacua 
with different a s the spin is 1/2 as before, but now the small- 
est change in the charge field is A(j) c = =bCi / 2 (tt/\^2), being 
C1/2 = I — 2 arccos[A/(2^3)]/7r the fractional charge of the 
excitation. As gs increases, Ci/ 2 decreases. This is consistent 
with a continuous transition to the MI phase, where the ele- 
mentary excitations are spinons without charge. In addition 
there is a pure charge excitation between the two possible val- 
ues of 4> c modulo 2tt for the same values of a c and a s . Its 
charge is Co = 1 — Ci/ 2 . 

B. The displacement operators 

The quantity z c L = (g\Ul\g) where [/£ = exp[i^X] 
where L is the length of the system and X = J2j J a ( n j] + 
riji)], was first proposed by Resta and Sorella as an indica- 
tor of localization in extended systems. 01 Using symmetry 
properties, Ortiz and one of us showed that for translation- 
ally invariant interacting systems with a rational number n/m 
of particles per unit cell the correct definition is (g\(Ul) m \g) 
13411 . This definition was used to characterize metal-insulator 
and metal- superconducting transitions in one dimensional lat- 
tice models. In the thermodynamic limit, L — > 00, |z£| is 
equal to 1 for a metallic periodic system and to for the 
insulating state. This is also true for non-interacting dis- 
ordered systems. O In the thermodynamic limit, z c L pro- 
vides a way to calculate the expectation value of the posi- 
tion operator (related to the macroscopic electric polarization 



b^ P = e(X)/L) and its fluctuations for a periodic system 
IS: 



lim 

l^oo znrn 

L 2 



-Im In z c r 



7c = P - P , 



lim 

L^oo (ZTTUl)' 



■lnl^l 2 = 



27T7TI 

(X 2 ) - {X)\ 



(20) 



where 7 C is the charge Berry phase E El H and 
Po some reference polarization. j c is defined modulo 2tt and 
therefore Po is defined modulo e/m. These results were ex- 
tended to z s L , where X is replaced by the difference between 
the position of up and down operators Y^j J a ( n j] ~ n jl)i P 
by the correspondin g di fference in polarization, and j c by the 
spin Berry phase ^. 1171 14411 For calculations in finite systems, 
more accurate results for j v are obtained calculating it directly 
rather than using the first Eq. (l20l> .l 171 13411 However, this ex- 
pression allows a trivial calculation of j c if all particles in the 
ground state are localized. fl7l1 In particular 7 C = for a CDW 
with maximum order parameter and j c = tt for a MI with one 
particle per site. By continuity and considering that j c can 
only be or ir modulo 2ir for a system with inversion sym- 
metry and one particle per site, these two values characterize 
the II and MI phases until the boundary of a phase transition 
is reached. 

The difficulty in finding a bosonized expression for U^ s is 
that the exponents are ill defined in a periodic system. Ac- 
tually this is the reason why P should be defined through 
the Berry phase or z£.0, [3 H3l Nevertheless, as noted 
earlier, 1 34] Z7£ is the operator that shifts all one particle mo- 
menta by -27r/L, to the next available wave vector to the 
left. Then, by inspection of the expressions of the fields Eqs. 
(1161171) . Z7£ acting on any state increases N- a by one and de- 
creases 7V +cr also by one. 14211 Consequently, Z7£ shifts U c (x) 
by —2^/2/L and leaves the other three fields invariant. From 
the commutation relation, 



it follows that 



[(j) c (x),Il c (y)} = iS(x - y), 



i / dx(f) c (x),U c (y)} 



(21) 



(22) 



Therefore the displacement operator ?7£ can be constructed as 
an exponential of the average charge field. A similar analysis 
can be carried out for U[ and both results can be written as 



Ul = exp(zv^), € = lf dx^(x). (23) 

In the region of parameters such that the interaction m v 
in the sine Gordon model Eq. (Tf4l is relevant according to 
the renormalization group, (j) v (x) gets frozen at the value that 
minimizes m v cos(VS(j) u ) and the calculation of J7£ becomes 
trivial. In particular, for gs± < (m c < corresponding to 
a CDW or II) we obtain $ a c = 7 C = 0, while in the MI phase 
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(U — > oo, m c > 0) the result is \fS<j>^ = 7c = tt, in agree- 
ment with the values anticipated above. iTTvIl The difference in 
polarization e/2 modulo e (see Eq. d20l ) is consistent with the 
transfer of half of the electrons in one lattice parameter which 
is required to convert a MI into a II in the strong coupling 
limit. 

In the BOI phase, there is a spontaneous breaking of the in- 
version symmetry in the thermodynamic limit and fractional 
values of 7 c /tt are allowed for A ^ 0. The analysis of 
the previous subsection indicates that in the BOI phase there 
are two non-equivalent values for the frozen charge field: 
V2(j)c = =b arccos[A/(2#3)] = ±7rCo/2. Then, using Eqs. 
(l20l and (l23l we obtain the following relation between the po- 
larization and the fractional charge of the elementary spinless 
excitation of the BOI phase: 

P-Po = ±^ c = ±^C . (24) 

The sign depends on the particular ground state out of the two- 
fold degenerated manifold which follows from the Z2 (spatial 
inversion) symmetry breaking. If P is chosen in such a way 
that P = for the MI (implying P = ±e/2), we also obtain: 

P = ±|c 1/2 . (25) 

V. RIGOROUS RESULTS ON AN APPROXIMATE 
EFFECTIVE HAMILTONIAN 

So far, no approximations have been made beyond that of 
the strong coupling limit. In this Section we describe a sim- 
plification of H e f f that allows to map it into an exactly solv- 
able model. We neglect the correlated hopping term t c h 
and treat the exchange term J as an independent parameter. 
Both approximations are usually applied to the t — J model 
after its derivation from the Hubbard model in the strong cou- 
pling limit Oli^] in spite of the fact that the correlated hop- 
ping term might be important for stabilizing a superconduct- 
ing resonance- valence-bond ground state. ll47H48ll 

To simplify the expression for the spin Hamiltonian and un- 
veil accidental symmetries that appear for particular combi- 
nations of the parameters, it is convenient to change the sign 
of the sf sf +1 and term in Hamiltonian © (this also 

changes the sign of the term (S^~ S~ +1 ) 2 -\-H.c. in Eq. (TTOl ) 
using the gauge transformation 

— > — for i = 4n and i = 4n + 1 
ctj — > -c] T for i = 4n + 2 and % — 4n + 1. (26) 

Then, neglecting t c h, adding the term in V in dTol . and tak- 
ing V = t, J = 2t, the spin 1 Hamiltonian takes the SU(2) 
invariant form: 

H eff = -*£[(S< • S l+1 ) 2 - 1] + £^T[(Sf) 2 - 1]. (27) 

i i 



For E = 0, this model has an SU(3) symmetry which is hid- 
den in this representation. The situation is somewhat similar 
to the supersymmetric t — J model for J = 2£.(4<| In 
our case, the SU(3) symmetry is made explicit usinga spin- 
particle transformation described in detail in Ref. 112 711 For 
E = 0, the model is equivalent to an isotropic SU(3) antifer- 
romagnetic Heisenberg model: 

H eff (E = 0)= tS^{j)S^{j + l\ (28) 

jeA,n,v 

where A is an arbitrary sublattice (J even or odd) in which 
the three relevant states transform under the fundamental rep- 
resentation S^ v ', while in the other sublattice B, the states 
transform under its conjugate representation S u ^. The Hamil- 
tonian Eq. (l28l is invariant under staggered conjugate SU(3) 
rotations, 1Z and TV, on sublattices A and B respectively and 
is integrable. 1I5TII52H53I1 The ground state of the exact solution 
is spin-dimerized and corresponds to the BOI in the original 
language. The degree of dimerization can be characterized by 
the order parameter 

£> = |^_ M -/i M +i|, (29) 

where hi-\^ = \i — 1, i — 1| is a projector on the SU(3) 
singlet spin state, \i — = \i — lm)\im) / y/3, at the 
bond (i — 1, i). Note that 3/^_i ? ^ coincides with the term with 
j = i — 1 in the isotropic SU(3) Heisenberg Hamiltonian Eq. 
(l28l . In the original fermion representation, \i — 1, i) has the 
form 

1* - 1. = + /ti T 4 - /tu4)i°)' ( 30 ) 

where k = i — [1 + (— 1)*]/2 is the site of energy — A/2 be- 
tween z and z + 1. A perfect dimerized state is constructed 
repeating this SU(3) singlet each two sites. There are two 
possibilities depending on the initial site chosen, reflecting 
the Z(2) symmetry breaking of the exact solution. The ex- 
act ground state has D = 0.42 16 Dq where Dq is the value of 
D for a perfect dimerized state. 1 5J] The value of the gap is 
rather small A = 0.173178t, and accordingly the correlation 
length £ = 21.0728505... is very large. B53I1 Therefore, one 
expects that very large systems should be studied to calculate 
numerically the bond-order parameter or other properties of 
the BOI. 

The significant degree of dimerization and the form of the 
SU(3) singlets Eq.(l30l» points out an important degree of co- 
valency in the BOI phase even in the strong coupling limit. 
This is in contrast to the solutions for the II and MI phases 
described previously for large V (see Section III). 

While the exact solution brings important insight into the 
physics of a dimerized state, one might wonder if the approx- 
imations made in the Hamiltonian, or the parameters chosen, 
render the results invalid for the IHM. If the IHM is described 
to first order in t, then t c h becomes irrelevant and the question 
is if the BOI survives after moving J and V from their values 
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at the exact solution to zero. The results presented in Sec- 
tions II and III, indicate that increasing V tends to close the 
BOI phase, while increasing J the effect is the opposite. The 
effect of small changes of J, V, and E from the exactly solv- 
able point on the energy of the dimer state, can be calculated 
by first order perturbation theory. In particular, because of the 
SU(3) symmetry of the ground state, the ground state expec- 
tation value of — (si ■ s^ + i — 77^+1/4), which is a projector 
over an SU(2) singlet with two components, is 2/3 times the 
expectation value of the SU(3) singlet with three components. 
We can write: 

SEboi = Eboi — E BOI = [- -(J — 2b) 

~\{V - t)](h iti+1 ) av + E((S? ) 2 - 1),(31) 

where E BOI is the energy at the exactly solvable point and 
(hi,i+i)av is the expectation value (ft^i+i) averaged over 
even and odd i. From the energy of the exact solution, we 
know that Q 

3(hi 9i+1 ) av = ((Si • S i+1 ) 2 - l) av = 1.796863..., (32) 

while because of symmetry (all spin 1 projections are equally 
probable) ((Sf) 2 ) = 2/3. 

Estimating the energy shifts for the other two phases as in 
Section II, neglecting terms of order t 2 we have 



with A] = (fl J+1] fl 3+n + 4- T 4 +U - /2V2W 
other ground state \d2) is obtained replacing the subscript 2j 
by 2j + 2 in the last two terms). M = L/(2a) is half the 
number of sites. Using the definitions in Section IV b and 
applying U L to our perfectly dimerized ground state \dl) we 
obtain: 

M-l 

n e fe ^(c*/t +lt /i +li +4 T /i +li -/i i 4 +1T )io>. 

(35) 

In the thermodynamic limit, (exp(i7r/M) + 2)/3 = 
exp(i7r/3M), the mean value of z L is: 

M-l 

z c L (dl) = (dl\U£\dl) ^ e i7T ] [ e i7T/3M = e*** . 

3=0 

This result implies a Berry phase equal to 47r/3. Assum- 
ing as before that the polarization of the MI phase is zero 
(Po = e/2 modulo e) we have for the polarization P = e/6, 
using Eq. (l20l> with m = 1. This result has a simple interpre- 
tation: the ideal dimer state is composed of isolated "diatomic 
molecules"and in this case P coincides with the result for one 
molecule. 1 34] In each molecule, 1/3 of the charge (the weight 
of the first term in each factor of Ea.(l34l)) has been displaced 
from site 2j to 2j + l, and since there is one molecule per two 
lattice parameters, the dipolar moment has increased by e/6. 



5E n = E u - E° n ~E-{V-t) 
5E M i = E M i -E° mi ~-(J- 2t) In 2. (33) 

Using these expressions to calculate the energies of the IHM 
(V = J = 0) for small E, one obtains that SEu = SEmi = 
1.39* when E/t = 1 - 2 In 2 ^ -0.39. For this value of E, 
Eq. (I3TI gives 5 Eboi — 1.1 3£ < SEu. Since we know that 
E BOI < E^j and E BOI < E® MI , this implies that the BOI 
phase is still that of lowest energy in the strong coupling limit 
of the IHM if E is decreased to ~ —0.3 9*. Numerical results 
suggest that while this shift in E is in the right direction, its 
magnitude is nearly two times smaller. JT2L I22L l23ll 



VI. POLARIZATION AND FRACTIONAL CHARGE 
EXCITATIONS 

In this section, we discuss the polarization and fractional 
charge excitations of the BOI phase. To render the discussion 
clearer, we begin by calculating these properties in a perfectly 
dimerized state with the above mentioned SU(3) symmetry. 
Due to the breaking of Z(2) symmetry, there are two such 
states as it is shown in Fig|Vl| Except for a trivial normal- 
ization constant, one of them can be written as: 

M-l 

\di)= n a >>, (34) 

j=0 




FIG. 1: Schematic plot of two dimerized ground states of the BOI. 

It is easy to generalize this result out of the SU(3) symmet- 
ric point and for the other perfectly dimerized state. The result 
is (modulo e) 

P(dl) = -P(d2) = |(m - 1) = |(1 - (^) 2 >. (36) 

The factors that multiply e/2 express the weight of doubly 
occupied sites, i.e., the sites with energy —A/2. Note that 
as U increases these factors decrease and they vanish when 
the system enters the MI phase. In the opposite limit, at the 
boundary with the II phase or inside it, there is also only one 
value of the polarization P = ±e/2 since it is defined modulo 
e. At this point, it is not possible to shift between two differ- 
ent polarized sates by application of an electric field and the 
system is not ferroelectric. Only inside the BOI phase one has 
ferroelectricity. 
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While we are unable to calculate the polarization in the 
ground state of the BOI phase, our previous result Eq. (l25l 
indicates that the result for the perfectly dimerized states Eq. 
(l36l holds in the general case. The charge of the elementary 
excitations of the BOI phase can be easily determined in the 
spin language. 1231 Fig. 2 shows a schematic representation of 
a spinon excitation in which the two possible ground states of 
the BOI phase are separated by a spin 1 at a given site i. This is 
not an eigenstate because the quantum fluctuations of the spin 
and the bond ordered regions are not included. However, the 
inclusion of these fluctuations does not change the charge of 
the soliton which is connecting both ground states. Since the 
original charge operator can be written as Q = -X~ 1) J \Sj I > 
if Sf = (representing either an empty site or a doubly oc- 
cupied one), the total change in spin is zero, and the charge 
of the excitation is Co = ((Sf) 2 ), where / is any site. If in- 
stead Sf = 1, the spin of the excitation is 1/2 and its charge is 
Ci/2 = (1 — (S[) 2 ). These results can be related with those 
of Section IV b and are consistent with Eqs. (l25l and (l36l . 



FIG. 2: Schematic plot of the soliton excitation of the dimerized 
ground state in the spin one representation (see Eq. i27l ). In the 
original fermionic representation, this excitation has spin S = 1/2 
and charge Ci/ 2 = 1/3 at the SU(3) symmetric point. 



spin-particle transformations introduced in Ref. 13 111 , we 
mapped the effective Hamiltonian into a spin one model. The 
spin language is useful to unveil hidden symmetries of the ef- 
fective model for particular combinations of the different pa- 
rameters. For instance, we showed that when the correlated 
hopping t c h is neglected and the exchange term J is treated as 
an independent parameter, the effective Hamiltonian becomes 
a biquadratic Heisenberg model with a single-ion anisotropy 
for V = t and J = 2t. In particular, the amplitude of the 
single-ion anisotropy vanishes in the region dominated by the 
charge transfer instability: U = A. At this point, the Hamil- 
tonian can be rewritten as an SU(3) antiferromagnetic Heisen- 
berg model using the general spin-particle transformations in- 
troduced in |22,[3(J]. The exact ground state of this model is a 
dimerized solution that corresponds to the BOI in the original 
fermionic language. In the spin one language, the low energy 
excitations of the dimerized solution are spinons. The S z = 
spinons carry zero spin and charge Cq = ±2/3 in terms of the 
original fermions. The S z = ±1 spinons have spin 5 = 1/2 
and charge Ci/ 2 = ±1/3 in the original fermionic version. 
Making U ^ A changes the value of Cq and C\ continuously 
as a function of U — A so the charge Q the solitons can also 
take irrational values. 

The above described results, which were obtained in the 
strong coupling limit, are in qualitative agreement with the 
weak coupling bosonization approach of Fabrizio et a/. lHHl 
The irrational values of the charge carried by each soliton are 
just a consequence of the asymmetry between the two sub- 
lattices introduced by the A term. 1 56] Since these excitations 
have a topological nature, it is expected that their characteris- 
tics will not depend on the interaction regime. 

In addition, using the bosonization procedure we showed 
that the polarization of the ferroelectric BOI is proportional 
to the charge Ci/ 2 of the elementary excitations (solitons). 
Therefore, the magnitude of the fractional charge carried by 
each soliton can be obtained experimentally by measuring the 
electric polarization of the ground state. 



VII. DISCUSSION 

In summary, we have derived an effective low energy 
Hamiltonian for the strong coupling regime of an extended 
IHM that includes a nearest-neighbor repulsion. Using the 
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